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This talk contains a short overview of the history of the interplay of the weak and 
the strong interaction and CP-violation. It describes the phenomenology and the 
basic physics mechanisms involved in the Standard Model calculations of K — > tttt 
decays with an emphasis on the evaluation of Penguin operator matrix-elements. 



1. Introduction 

In this conference in honour of Arkady Vainshtein's 60th birthday, a dis- 
cussion of the present state of the art in analytical calculations of relevance 
for Kaon decays is very appropriate given Arkady's large contributions to 
the field. He has summarised his own contributions on the occasion of ac- 
cepting the Sakurai Prize. 1 This also contains the story of how Penguins, 
at least the diagram variety, got their name. In Fig. 1 I show what a real 
(Linux) Penguin looks like and the diagram in a "Pcnguinized" version. 
This talk could easily have had other titles, examples are "QCD and Weak 
Interactions of Light Quarks" or 'Penguins and Other Graphs." In fact 
I have left out many manifestations of Penguin diagrams. In particular I 
do not cover the importance in B decays where Penguins were first exper- 
imentally verified via B — > K^j, but are at present more considered a 
nuisance and often referred to as "Penguin Pollution." Penguins also play 
a major role in other Kaon decays, reviews of rare decays where pointers 
to the literature can be found are Refs [2,3,4]. 

In Sect. 2 I give a very short historical overview. Sects 3 and 4 discuss 
the main physics issues and present the relevant phenomenology of the 
AI =1/2 rule and the CP-violating quantities e and e'/e. The underlying 
Standard Model diagrams responsible for CP-violation are shown there as 
well. The more challenging part is to actually evaluate these diagrams in 
the presence of the strong interaction. We can distinguish several regimes of 
momenta which have to be treated using different methods. An overview 



1 



2 




Figure 1. The comparison between a (Linux) Penguin and the Penguin diagram. (Linux 
Penguin from Neal Tucker, (http://www.isc.tamu.edu/~lewing/linux/).) 

is given in Sect. 5 where also the short-distance part is discussed. The 
more difficult long-distance part has a long history and some approaches 
are mentioned, but only my favourite method, the A-boson or fictitious 
Higgs exchange, is described in more detail in Sect. 6 where I also present 
results for the main quantities. For two particular matrix-elements, those 
of the electroweak Penguins Q-j and Q%, a dispersive analysis allows to 
evaluate these in the chiral limit from experimental data. This is discussed 
in Sect. 7. We summarise our results in the conclusions and compare with 
the original hopes from Arkady and his collaborators. This talk is to a large 
extent a shorter version of the review [5]. 

2. A short historical overview 

The weak interaction was discovered in 1896 by Becquerel when he discov- 
ered spontaneous radioactivity. The next step towards a more fundamental 
study of the weak interaction was taken in the 1930s when the neutron was 
discovered and its /3-decay studied in detail. The fact that the proton and 
electron energies did not add up to the total energy corresponding to the 
mass of the neutron, made Pauli suggest the neutrino as a solution. Fermi 
then incorporated it in the first full fledged theory of the weak interaction, 
the famous Fermi four-fermion 6 interaction. 
Gf 

£ Fermi = ~7S (1 ~ 7s) n] [ej^ (1 - 75) v\ . (1) 
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The first fully nonhadronic weak interaction came after world-war two 
with the muon discovery and the study of its ,5-decay. The analogous 
Lagrangian to Eq. (1) was soon written down. At that point T.D. Lee and 
C.N. Yang 7 realized that there was no evidence that parity was conserved 
in the weak interaction. This quickly led to a search for parity violation 
both in nuclear decays 8 and in the decay chain tt + — > — > e+^P^i/^. 9 
Parity violation was duly observed in both cases. These experiments and 
others led to the final form of the Fermi Lagrangian given in Eq. (1). 10,11 

During the 1950s steadily more particles were discovered providing many 
puzzles. These were solved by the introduction of strangeness, 12 ' 13 of what 
is now known as the Kl and the K$ 14 and the "eightfold way" of classifying 
the hadrons into symmetry- multiplets. 15 

Subsequently Cabibbo realized that the weak interactions of the strange 
particles were very similar to those of the nonstrange particles. 16 He pro- 
posed that the weak interactions of hadrons occurred through a current 
which was a mixture of the strange and non-strange currents with a mixing 
angle now universally known as the Cabibbo angle. The hadron symmetry 
group led to the introduction of quarks 17 as a means of organising which 
SU(3)v multiplets were present in the spectrum. 

In the same time period the Kaons provided another surprise. Mea- 
surements at Brookhaven 18 indicated that the long-lived state, the Kl, 
did occasionally decay to two pions in the final state as well, showing that 
CP was violated. Since the CP-violation was small, explanations could be 
sought at many scales, an early phenomenological analysis can be found in 
Rcf. [19], but as the socalled superweak model 20 showed, the scale of the 
interaction involved in CP- violation could be much higher. 

The standard model for the weak and electromagnetic interactions of 
lcptons was introduced in the same period. The Fermi theory is nonrenor- 
malisable. Alternatives based on Yang-Mills 21 theories had been proposed 
by Glashow 22 but struggled with the problem of having massless gauge 
bosons. This was solved by the introduction of the Higgs mechanism by 
Weinberg and Salam. The model could be extended to include the weak 
interactions of hadrons by adding quarks in doublets, similar to the way the 
lcptons were included. One problem this produced was that loop-diagrams 
provided a much too high probability for the decay Kl — ► com- 
pared to the experimental limits. These socalled flavour changing neutral 
currents (FCNC) needed to be suppressed. The solution was found in the 
Glashow-Iliopoulos-Maiani mechanism. 23 A fourth quark, the charm quark, 
was introduced beyond the up, down and strange quarks. If all the quark 
masses were equal, the dangerous loop contributions to FCNC processes 
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cancel, the socalled GIM mechanism. This allowed a prediction of the 
charm quark mass, 24 soon confirmed with the discovery of the J/ip- 

In the mean time, QCD was formulated. 25 The property of asymptotic 
freedom 26 was established which explained why quarks at short distances 
could behave as free particles and at the same time at large distances be 
confined inside hadrons. 

The study of Kaon decays still went on, and an already old problem, 
the AI = 1/2 rule saw the first signs of a solution. It was shown 27,28 
that the short-distance QCD part of the nonlcptonic weak decays provided 
already an enhancement of the AI = 1/2 weak AS — 1 transition over the 
AI = 3/2 one. The ITEP group extended first the Gaillard-Lee analysis 
for the charm mass, 29 but then realized that in addition to the effects 
that were included in Refs [27,28], there was a new class of diagrams that 
only contributed to the AI =1/2 transition. 30,31 While, as we will discuss 
in more detail later, the general class of these contributions, the socalled 
Penguin-diagrams, is the most likely main cause of the AI = 1/2 rule, the 
short-distance part of them provide only a small enhancement contrary to 
the original hope. A description of the early history of Penguin diagrams, 
including the origin of the name, can be found in the 1999 Sakurai Prize 
lecture of Vainshtein. 1 

Penguin diagrams at short distances provide nevertheless a large amount 
of physics. The origin of CP-violation was (and partly is) still a mystery. 
The superweak model explained it, but introduced new physics that had no 
other predictions. Kobayashi and Maskawa 32 realized that the framework 
established by Ref. [23] could be extended to three generations. The really 
new aspect this brings in is that CP-violation could easily be produced 
at the weak scale and not at the much higher superweak scale. In this 
Cabibbo-Kobayashi-Maskawa (CKM) scenario, CP-violation comes from 
the mixed quark-Higgs sector, the Yukawa sector, and is linked with the 
masses and mixings of the quarks. Other mechanisms at the weak scale 
also exist, as e.g. an extended Higgs sector. 33 

The inclusion of the CKM mechanism into the calculations for weak 
decays was done by Oilman and Wise 34 ' 35 which provided the prediction 
that e'/e should be nonzero and of the order of 10~ 3 . Guberina and Pec- 
cei 36 confirmed this. This prediction spurred on the experimentalists and 
after two generations of major experiments, NA48 at CERN and KTeV at 
Fermilab have now determined this quantity and the qualitative prediction 
that CP-violation at the weak scale exists is now confirmed. Much stricter 
tests of this picture will happen at other Kaon experiments as well as in B 
meson studies. 
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The K°-K a mixing has QCD corrections and CP-violating contribu- 
tions as well. The calculations of these required a proper treatment of 
box diagrams and inclusions of the effects of the AS = 1 interaction 
squared. This was accomplished at one-loop by Gilman and Wise a few 
years later. 37 ' 38 

That Penguins had more surprises in store was shown some years later 
when it was realized that the enhancement originally expected on chiral 
grounds for the Penguin diagrams 30 > 31 was present, not for the Penguin 
diagrams with gluonic intermediate states, but for those with a photon. 39 
This contribution was also enhanced in its effects by the AI = 1/2 rule. 
This lowered the expectation for e'/e, but it became significant after it was 
found that the top quark had a very large mass. Flynn and Randall 40 
reanalysed the electromagnetic Penguin with a large top quark mass and 
included also Z° exchange. The final effect was that the now rebaptized 
electroweak Penguins could have a very large contribution that could even 
cancel the contribution to e'/e from gluonic Penguins. This story still 
continues at present and the cancellation, though not complete, is one of 
the major impediments to accurate theoretical predictions of e'/e. 

The first calculation of two-loop effects in the short-distance part was 
done in Rome 41 in 1981. The value of Aq^p has risen from values of 
about 100 MeV to more than 300 MeV. A full calculation of all operators 
at two loops thus became necessary, taking into account all complexities 
of higher order QCD. This program was finally accomplished by two inde- 
pendent groups. One in Munich around A. Buras and one in Rome around 
G. Martinelli. 



3. K — > 7T7T and the AI — 1/2 rule 

The underlying qualitative difference we want to understand is the AI = 
1/2 rule. We can try to calculate K — > tttt decays by simple W + exchange. 
For K + — ► 7T + 7r° we can draw the two Feynman diagrams of Fig. 2(a). The 
VF + -hadron couplings are known from semi-leptonic decays. This approxi- 
mation agrees with the measured decay within a factor of two. 

A much worse result appears when we try the same for K° -» 7r°7r°. As 
shown in Fig. 2(b) there is no possibility to draw diagrams similar to those 
in Fig. 2(a). The needed vertices always violate charge-conservation. So 
we expect that the neutral decay should be small compared with the ones 
with charged pions. Well, if we look at the experimental results, we see 

T(K a — > ir°ir a ) = \?{Ks — > ^°7r°) = 2.3 • 10~ 12 MeV 
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K° 



(a) 




W+ 



(b) 



Figure 2. (a) The two naive VK+-exchange diagrams for A"+ > 7r+7r . (b) No simple 

W+-exchange diagram is possible for K° > 7r 7r°. 



T(K + — > tt+tt") = 1.1 • 1(T 14 MeV (2) 

So the expected zero one is by far the largest !!! 

The same conundrum can be expressed in terms of the isospin ampli- 
tudes: a 

A[K° -> vr 7r ] = y/T/3Ao - ^/2/3A 2 
A[K° -> tx+t:-] = yAjSAa + ^/TjQA 2 
A[K+ -» tt+tt ] = (V3/2)A 2 . (3) 

The above quoted experimental results can now be rewritten as 

IA>M 2 | exp = 22 (4) 

while the naive M /+ -exchange discussed would give 

l^oM 2 | naive = V2. (5) 

This discrepancy is known as the problem of the AI =1/2 rule. 

Some enhancement comes from final state 7T7r-rescattering. Removing 
these and higher order effects in the light quark masses one obtains 42 ' 43 

\A /A 2 \ x = l7.8. (6) 

This changes the discrepancy somewhat but is still different by an order 
of magnitude from the naive result (5). The difference will have to be 
explained by pure strong interaction effects and it is a qualitative change, 
not just a quantitative one. 

We also use amplitudes without the final state interaction phase: 

Ai = -ia ie l&1 (7) 



a The sign convention is the one used in the work by J. Prades and myself. 
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for I = 0, 2. 5j is the angular momentum zero, isospin I scattering phase 
at the Kaon mass. 



4. K — > 7T7T, e, e' /e 

The K°, K° states have sd, ds quark content. CP acts on these states as 

CP\K°) = -\K°) . (8) 
We can construct eigenstates with a definite CP transformation: 



K° 



j= (k° (+)K°) , CP\K m = +(-)\K m . (9) 



The main decay mode of K°-\ike states is tttt. A two pion state with charge 
zero in spin zero is always CP even. Therefore the decay K\ — * tttt is 
possible but K 2 — » nn is impossible: K 2 — » 7T7r7r is possible. Phase-space for 
the 7T7T decay is much larger than for the three-pion final state. Therefore 
if we start out with a pure K° or K° state, the K 2 component in its wave- 
function lives much longer than the Ki component such that after a long 
time only the K 2 component survives. 

In the early sixties, as you see it pays off to do precise experiments, one 
actually measured 18 

showing that CP is violated. This leaves us with the questions: 

??? Does K\ turn in to K 2 (mixing or indirect CP violation)? 
??? Does K 2 decay directly into tttt (direct CP violation)? 

In fact, the answer to both is YES and is major qualitative test of the stan- 
dard model Higgs-fermion sector and the CftTM-picture of CP- violation. 

The presence of CP- violation means that K\ and K 2 are not the mass 
eigenstates, these are 

K s(L) = ^==(K 1(2) +eK 2(1) ) . (11) 

They are not orthogonal since the Hamiltonian is not hermitian. 
We define the observables 

A(K L -> (7T7T) J=0 ) 

£ ~ A(K S -» (7r7r) 7=0 ) 

1 ( A(K L ^(™) I=2 ) A(K S ^ (™) I=2 ) \ 
V2\A(K s ^(nn) I=0 ) A(K S - (^)/ =0 ) J ' 1 ' 
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( a ) (b) 

Figure 3. (a) The box diagram contribution to K°K° mixing. Crossed versions and 
diagrams with extra gluons etc. are not shown, (b) The Penguin diagram contribution 
to K — > ■kit. Extra gluons and crossed versions etc. are not shown. 



The latter has been specifically constructed to remove the K°-K° transi- 
tion. |e| is a directly measurable as ratios of decay rates. 

We now make a series of experimentally valid approximations, 

|Ima |,|Ima 2 | « |Rea 2 | « |Rea |, |e|,|e|<<l, \e'\ « \e\ , (13) 

to obtain the usually quoted expression 

£ , = J_ e i(5 2 -5 D ) R£^2 / Ima2 _ Imao \ _ (u) 
y/2 Reao \ Rea2 Reao / 

Exp er iment ally, 44 

|e| = (2.271 ± 0.017) • 10~ 3 . (15) 

The set of diagrams, depicted schematically in Fig. 3(a), responsible for 
K°K° mixing are known as box diagrams. It is the presence of the virtual 
intermediate quark lines of up, charm and top quarks that produces the 
CP-violation. 

The experimental situation on e' /e was unclear for a long time. Two 
large experiments, NA31 at CERN and E731 at FNAL, obtained conflict- 
ing results in the mid 1980's. Both groups have since gone on and build 
improved versions of their detectors, NA48 at CERN and KTeV at FNAL. 
e'/e is measured via the double ratio 

Ke U J 6 l 1 T(K L -> n\°)/T(K s - ttV) / ' 

The two main experiments follow a somewhat different strategy in measur- 
ing this double ratio, mainly in the way the relative normalisation of Kl 
and Kg components is treated. After some initial disagreement with the 
first results, KTeV has reanalysed their systematic errors and the situation 
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Table 1. Recent results on e'/e. The 
years refer to the data sets. 



NA31 


(23.0 ± 6.5) 


X 


io- 


-4 


E731 


(7.4 ±5.9) 


X 


10" 


4 


KTeV 96 


(23.2 ±4.4) 


X 


io- 


-4 


KTcV 97 


(19.8 ±2.9) 


X 


io- 


-4 


NA48 97 


(18.5 ± 7.3) 


X 


io- 


-4 


NA48 98+99 


(15.0 ±2.7) 


X 


io- 


-4 


ALL 


(17.2 ±1.8) 


X 


10" 


-4 



for e'/e is now quite clear. We show the recent results in Table 1. The 
data are taken from Ref. [45] and the recent reviews in the Lepton-Photon 
conference. 46,47 

The Penguin diagram shown in Fig. 3(b) contributes to the direct CP- 
violation as given by e'. Again, py-couplings to all three generations show 
up so CP-violation is possible in K — > nn. This is a qualitative prediction 
of the standard model and borne out by experiment. The main problem 
is now to embed these diagrams and the simple FT-exchange in the full 
strong interaction. The AI = 1/2 rule shows that there will have to be 
large corrections to the naive picture. 

5. From Quarks to Mesons: a Chain of Effective Field 
Theories 

The full calculation in the presence of the strong interaction is quite dif- 
ficult. Even at short distances, due to the presence of logarithms of large 
ratios of scales, a simple one-loop calculation gives very large effects. These 
need to be resummed which fortunately can be done using renormalisation 
group methods. 

The three steps of the full calculation are depicted in Fig. 4. First we 
integrated out the heaviest particles step by step using Operator Product 
Expansion methods. The steps OPE we describe in the next subsections 
while step ??? we will split up in more subparts later. 

5.1. Step I: from SM to OPE 

The first step concerns the standard model diagrams of Fig. 5(a). We 
replace their effect with a contribution of an effective Hamiltonian given by 

W eff - E CiWQiin) = ^V ud V: s E [zi - Vi^r^t) Qi ■ (17) 

In the last part we have real coefficients z,- L and and the CKM-matrix- 
elements occurring are shown explicitly. The four-quark operators Q% can 
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ENERGY SCALE 



FIELDS 



Effective Theory 



Mi 



< 



M 



K 



W,Z,j,g; 
t, b, c, s, u, d 



J| using OPE 



7,5; (i,e,i>e; 
s, d, u 



7; M,e,^; 

7T, K, 1] 



Standard Model 



QCD,QED, 

t,|AS| = 1,2 



CHPT 



Figure 4. A schematic exposition of the various steps in the calculation of nonleptonic 
matrix-elements. 





(b) 



(a) 

Figure 5. (a) The standard model diagrams to be calculated at a high scale, (b) The 
diagrams needed for the matrix-elements calculated at a scale /ij/ f» m-w using the 
effective Hamiltonian. 



be found in e.g. Ref. [48]. 

We calculate now matrix-elements between quarks and gluons in the 
standard model using the diagrams of Fig. 5(a) and equate those to the 
same matrix-elements calculated using the effective Hamiltonian of Eq. (17) 



11 



Table 2. The Wilson coefficients and their main source at the scale fin = mw m the 
NDR-scheme. 



21 


0.053 


9,7-h>ox 


2/6 


-0.0019 


g-Pcnguin 


22 


0.981 


W + -cxchange g, 7-box 


2/7 


0.0009 


7, Z-Penguin 


2/3 


0.0014 


g, Z-Penguin WW-box 


2/8 


0. 




2/4 


-0.0019 


g-Pcnguin 


2/9 


-0.0074 


7, Z-Penguin Wiy-box 


2/5 


0.0006 


g-Pcnguin 


2/10 


0. 





and the diagrams of Fig. 5(b). This determines the value of the Zi and yi. 
The top quark and the W and Z bosons are integrated out all at the same 
time. There should be no large logarithms present due to that. The scale 
ji = [in in the diagrams of Fig. 5(b) of the OPE expansion diagrams should 
be chosen of the order of the W mass. The scale nw in the Standard Model 
diagrams of Fig. 5(a) should be chosen of the same order. 
Notes : 

• In the Penguin diagrams CP-violation shows up since all 3 generations 
are present. 

• The equivalence is done by calculating matrix-elements between Quarks 
and Gluons 

• The SM part is /^-independent to a s (nw)- 

• OPE part: The dependence of C^h) cancels the fin dependence of 
the diagrams to order a| (/!#). 

This procedure gives at nw = I^h = Mw in the NDR-scheme b the 
numerical values given in Table 2. In the same table I have given the main 
source of these numbers. Pure tree-level VF-exchange would have only given 
Z2 = 1 and all others zero. Note that the coefficients from 7, Z exchange 
are similar to the gluon exchange ones since as at this scale is not very big. 



5.2. Step II 

Now comes the main advantage of the OPE formalism. Using the renormal- 
isation group equations we can calculate the change with \x of the Cj, thus 
rcsumming the log (m^//x 2 ) effects. The renormalisation group equations 
(RGEs) for the strong coupling and the Wilson coefficients are 

M^5s(m) = /?(5s(m)). ^-^C i (fj,)=-f ji (gs(tJ,),a)C j (fj,). (18) 

(3 is the QCD beta function for the running coupling. The coefficients 7^ 
are the elements of the anomalous dimension matrix 7. They can be derived 

b The precise definition of the four-quark operators Qi comes in here as well. See the 
lectures by Buras 49 for a more extensive description of that. 
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Table 3. The Wilson coefficients Zi and yt at a scale AtQPE = ^00 MeV 
in the NDR scheme and in the X-boson scheme at fix = 900 MeV. 



i 




Vi 


Zi 


Vi 




HOPE = 


0.9 GeV 


MX = 


0.9 GeV 


21 


-0.490 


0. 


-0.788 


0. 


22 


1.266 


0. 


1.457 


0. 


23 


0.0092 


0.0287 


0.0086 


0.0399 


24 


-0.0265 


-0.0532 


-0.0101 


-0.0572 


25 


0.0065 


0.0018 


0.0029 


0.0112 


26 


-0.0270 


-0.0995 


-0.0149 


-0.1223 


27 


2.6 10~ 5 


-0.9 10~ 5 


0.0002 


-0.00016 


28 


5.3 10~ 5 


0.0013 


6.8 10~ 5 


0.0018 


29 


5.3 10~ 5 


-0.0105 


0.0003 


-0.0121 


210 


-3.6 10~ 5 


0.0041 


-8.7 10~ 5 


0.0065 



from the infinite parts of loop diagrams and this has been done to one 50 
and two loops. 51 The series in a and as is known to 



, n as ,i / as\ 2 „ ot „ aa a ,_„. 
^if+Mif) + - (19) 

Many subtleties are involved in this calculation. 49,51 They all are related to 
the fact that everything at higher loop orders need to be specified correctly, 
and many things which are equal at tree level are no longer so in d ^ 4 and 
at higher loops, see the lectures [49] or the review [52]. The numbers below 
are obtained by numerically integrating Eq. (18). 53,54 

We perform the following steps to get down to a scale Hope around 
1 GeV. Starting from the Zi and j/j at the scale fin : 

(1) solve Eqs. (18); run from fin to fib- 

(2) At fj,b~ rrib remove 6-quark and match to the theory without b by 
calculating matrix-elements of the effective Hamiltonian in the five and in 
the four-quark picture and putting them equal. 

(3) Run from fib to fi c K m c . 

(4) At fi c remove the c-quark and match to the theory without c. 

(5) Run from fi c to /^OPE- 

Then all large logarithms including mw, m z, "it, mb and m c , are summed. 

With the inputs m t {m t ) = 166 GeV, a = 1/137.0, a s (m z ) = 0.1186 
which led to the initial conditions shown in Table 2, we can perform the 
above procedure down to hope- Results for 900 MeV are shown in columns 
two and three of Table 3. z\ and Z2 have changed much from and 1. This 
is the short-distance contribution to the AI =1/2 rule. We also see a large 
enhancement of and yg, which will lead to our value of e'. 
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5.3. Step III: Matrix- elements 

Now remember that the Ci depend on [ioPE (scale dependence) and on 
the definition of the Qi (scheme dependence) and the numerical change 
in the coefficients due to the various choices for the Qi possible is not 
negligible. It is therefore important both from the phenomenological and 
fundamental point of view that this dependence is correctly accounted for 
in the evaluation of the matrix-elements. We can solve this in various ways. 

• Stay in QCD Lattice calculations. 55 

• ITEP Sum Rules or QCD sum rules. 56 

• Give up =^> Naive factorisation. 

• Improved factorisation 

• X-boson method (or fictitious Higgs method) 

• Large N c (in combination with something like the X-boson method.) 
Here the difference is mainly in the treatment of the low-energy hadronic 
physics. Three main approaches exist of increasing sophistication. 

© CHPT: As originally proposed by Bardeen-Buras-Gerard 57 and now 
pursued mainly by Hambye and collaborators. 58 

© ENJL (or extended Nambu-Jona-Lasinio model 59 ): As mainly done 
by myself and J. p ra des. 60 > 48 > 61 > 53 > 54 

© LMD or lowest meson dominance approach. 62 These papers stay with 
dimensional regularisation throughout. The A-boson corrections dis- 
cussed below, show up here as part of the QCD corrections. 

• Dispersive methods Some matrix-elements can in principle be deduced 
from experimental spectral functions. 

Notice that there other approaches as well, e.g. the chiral quark model. 63 
These have no underlying arguments why the ^-dependence should cancel, 
but the importance of several effects was first discussed in this context. I 
will also not treat the calculations done using bag models and potential 
models which similarly do not address the /x-dependence issue. 

6. The X-boson Method and Results using ENJL for the 
Long Distance 

We want to have a consistent calculational scheme that takes the scale and 
scheme dependence into account correctly. Let us therefore have a closer 
look at how we calculate the matrix-elements using naive factorisation. We 
start from the four-quark operator: 



c Which of course means that calculations exist only for simpler matrix-elements for the 
more sophisticated approaches. 
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Xb _ I tsO 



o 
o 

K_ {> \ gluono 
o 
o 



large q x 



(a) (b) 

Figure 6. (a) The leading in 1/N C contribution from Xb exchange, (b) The large 
momentum part of the Xb exchange matrix-element. 

=> See it as a product of currents or densities. 

=> Evaluate current matrix-elements in low energy theory or model or from 
experiment. 

=> Neglect extra momentum transfer between the current matrix elements. 
The main lesson here is that currents and densities are easier to deal with. 
We also need to go beyond the approximation in the last step. To obtain 
well defined currents, we replace the four-quark operators by exchanges of 
fictitious massive X-bosons coupling to two-quark currents or densities. 

H = Y J Ci(iiOP E )Qi => Y,9 iXiJi - ( 2 °) 

i i 

• This is a well defined scheme of nonlocal operators. 

• The matching to obtain the coupling constants gi from the Cj is done 
with matrix-elements of quarks and gluons. 

A simple example is the one needed for the Bk parameter. The four- 
quark operator is replaced by the exchange of one X-boson Xb'- 

C(m)^7m( 1 - Ibjsd-i^l - 7 5 )s =*> gBX^dj^l - j 5 )s . (21) 

Taking a matrix-element between quarks at next-to-leading order in as 
gives 

C(hope) (1 + ras(nop E )) - (<?I/M|J (l + r'a s + r" log (Ml B / M 2 )) .(22) 

The coefficients r and r' take care of the scheme dependence. The l.h.s. is 
scale independent to the required order in aj. The effect of these coeffi- 
cients surprisingly always went in the direction to improve agreement with 
experiment 48 > 54 as can be seen from columns 4 and 5 in Table 3. 

The final step is the matrix-element of Xs-boson exchange. For this we 
split the integral over the Xb momentum qx in two parts 

dq\ . (23) 



fdq x = 




I 




Jo 


J n 
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The leading in iV c contribution is depicted in Fig. 6(a) and corresponds to 
the large N c factorisation. The large momentum regime is evaluated by the 
diagram in Fig. 6(b), since the large momentum must flow back through 
quarks and gluons. Hadronic exchanges are power suppressed because of 
the form factors involved. The as present already suppresses by iV c so the 
matrix-element of this part can be evaluated using factorisation. This part 
cancels the r" \og(M XB / y?) present in (22). The final part with small qx 
momentum in the integral then needs to be evaluated nonperturbatively. 
Here one can use Chiral Perturbation Theory, the ENJL model or meson 
exchange approximations with various short-distance constraints. 

Let me now show some results from Refs [48,54]. The chiral limit cou- 
pling G*8 responsible for the octet contribution, it is 1 in the naive approx- 
imation and about 6 when fitted to experiment, 42,43 is shown in Fig. 7. As 
can be seen, the matching between the short and long distance is reasonable 
both for the real and imaginary parts. The value of ReGs is dominated by 
the matrix-element of Q\ and Qi but about 30-60% comes from the long 
distance Penguin part of Q2- The result for the ImGs corresponds to a 
value of the Qe matrix-element much larger than usually assumed. We 
obtained Bq « 2-2.5 while it is usually assumed to be less than 1.5. 

Putting our results in (14) we obtain a chiral limit value for (e'/e) x of 
about 6 • 10 -3 . 

We now add the main isospin breaking component 64 and the effect 
of final state interaction (FSI). 65 The latter in our case has mainly effect on 
the forefactor Rec^/Reao in (14) since the ratios of imaginary parts have 
been evaluated to the same order in p 2 in CHPT and thus receive no FSI 
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Table 4. The values of the VEVs in the NDR scheme at fi R = 2 GeV. The most 
recent dispersive results are line 3 to 5. The other results are shown for comparison. 
Errors are those quoted in the papers. Adapted from Ref. [71]. 



Reference 


{0\Og>\0)» D * 




B 7 = Bg = 1 


-(5.4 ±2.2) ■ 10~ 5 GeV 6 


(1.0 ± 0.4) ■ 10- 3 GeV 6 


Bijnens et al. [71] 
Knccht ct al. [68] 
Cirigliano et al. [70] 


-(4.0 ±0.5) ■ 10" 5 GeV 6 
-(1.9 ±0.6) ■ 10~ 5 GeV 6 
-(2.7 ± 1.7) ■ 10" 5 GeV 6 


(1.2 ± 0.5) ■ 10" 3 GeV 6 
(2.3 ± 0.7) ■ 10~ 3 GeV 6 
(2.2 ± 0.7) ■ 10" 3 GeV 6 


Donoghuc ct al.[67] 
Narison [69] 
lattice [72] 
ENJL [54] 


-(4.3 ±0.9) ■ W~ b GeV 6 
-(3.5 ±1.0) • 10" 5 GeV 6 
-(2.6 ±0.7) ■ 10~ 5 GeV 6 
-(4.3 ±0.5) ■ 10~ 5 GeV 6 


(1.5 ± 0.4) ■ 10" 3 GeV 6 
(1.5 ± 0.3) ■ 10~ 3 GeV 6 
(0.74 ±0.15) ■ 10~ 3 GeV 6 
(1.3 ± 0.2) ■ 10~ 3 GeV 6 



corrections. The final result is 54 e'/s w 1.5 ■ 10 3 with an error ^ 50%. 
7. Dispersive Estimates for {Q7) and (Qs) 

Some of the matrix-elements we want can be extracted from experimental 
information in a different way. The canonical example is the mass difference 
between the charged and the neutral pion in the chiral limit which can be 
extracted from a dispersive integral over the difference of the vector and 
axial vector spectral functions. 66 

This idea has been pursued in the context of weak decay in a scries 
of papers by Donoghue, Golowich and collaborators. 67 The matrix-element 
of Q-j could be extracted directly from these data. To get at the matrix- 
element of Qs is somewhat more difficult. Ref. [67] extracted it first by 
requiring /i-indcpendcncc, this corresponds to extracting the matrix ele- 
ment of Qs from the spectral functions via the coefficient of the dimension 
6 term in the operator product expansion of the underlying Green's func- 
tion. The most recent papers using this method are Refs. [68,69,70] and [ 
71]. In the last two papers also some QCD corrections were included which 
had a substantial impact on the numerical results. 

The results are given in Table. 4. The operator is related by a 

(2) 

chiral transformation to Q7 and 6 to Qs. The numbers are valid in 
the chiral limit. The various results for the matrix-clement of Og arc in 
reasonable agreement with each other. The underlying spectral integral, 
evaluated directly from data in Refs. [69], [70] and [71], or via the minimal 
hadronic ansatz [68] arc in better agreement. The largest source of the 
differences is the way the different results for the underlying evaluation of 
Og 2 ' come back into O^. 

The results for are also in reasonable agreement. Ref. [71] uses two 
approaches. First, the matrix-element for O e can be extracted via a sim- 
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ilar dispersive integral over the scalar and pseudoscalar spectral functions. 
The requirements of short-distance matching for this spectral function com- 
bined with a saturation with a few states imposes that the nonfactorisable 
part is suppressed and the number and error quoted follows from this. Ex- 
tracting the coefficient of the dimension 6 operator in the expansion of the 
vector and axial-vector spectral functions yields a result comparable but 
with a larger error of about 0.9. Ref. [68] uses a derivation based on a sin- 
gle resonance plus continuum ansatz for the spectral functions and assumes 
a typical large N c error of 30%. This ansatz worked well for lower moments 
of the spectral functions which can be tested experimentally. Adding more 
resonances allows for a broader range of results. Ref. [70] chose to enforce 
all the known constraints on the vector and axial-vector spectral functions 
to obtain a result. This resulted in rather large cancellations between the 
various contributions making an error analysis more difficult. A reasonable 
estimate lead to the value quoted. 

The reason why the central value based on the same data can be so 
different is that the quantity in question is sensitive to the energy regime 
above 1.3 GeV where the accuracy of the data is rather low. 

8. Conclusions 

Penguins are alive and well, they provide a sizable part of the AI =1/2 en- 
hancement though mainly through long distance Penguin like topologies in 
the evaluation of the matrix- element of Q2- They have found a much richer 
use in the CP violation phenomenology. For the electroweak Penguins, 
calculations are in qualitative agreement but more work is still needed to 
get the errors down. For the strong Penguins, the work I have presented 
here shows a strong enhancement over factorisation with Bq significantly 
larger than one. The latter conclusion is similar to the one derived from the 
older more phenomenological arguments where the coefficients were taken 
at a low scale and the matrix-elements for Qq taken from the the value of 
the AI = 1/2 rule. This also indicated a rather large enhancement of the 
matrix-element of Qq over the naive factorisation. 
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